Noise Induced Intermittency in a Superconducting Microwave Resonator 
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We experimentally and numerically study a NbN superconducting stripline resonator integrated 
with a microbridge. We find that the response of the system to monochromatic excitation exhibits 
intermittency, namely, noise-induced jumping between coexisting steady-state and limit-cycle re- 
sponses. A theoretical model that assumes piecewise linear dynamics yields partial agreement with 
the experimental findings. 



Nonlinear response of superconducting RF devices can 
be exploited for a variety of applications such as noise 
squeezing [l[ , bifurcation amplification 0, 0] and res- 
onant readout of qubits 0]. Recently we have reported 
on an instability found in NbN superconducting stripline 
resonators in which a short section of the stripline was 
made relatively narrow, forming thus a microbridge @,[3|- 
In these experiments a monochromatic pump tone, hav- 
ing a frequency close to one of the resonance frequencies, 
is injected into the resonator and the reflected power off 
the resonator is measured. We have discovered that there 
is a certain zone in the pump frequency - pump ampli- 
tude plane, in which the resonator exhibits limit-cycle 
(LC) response resulting in self-sustained modulation of 
the reflected power. Moreover, to account for the ex- 
perimental findings we have proposed a simple piecewise 
linear model, which attributes the resonator's nonlinear 
response to thermal instability occurring in the micro- 
bridge [8|. In spite of its simplicity, this model yields 
a rich variety of dynamical effects. In particular, as we 
show below, it predicts the occurrence of intermittency, 
namely the coexistence of different LC and steady state 
solutions, and noise-induced jumping between them. In 
the present paper we study both theoretically and ex- 
perimentally noise-induced transitions between different 
metastable responses. We employ a ID map to find the 
possible LC solutions of the system and to find conditions 
for the occurrence of intermittency. Experimentally we 
present measurements showing both, intermittency be- 
tween a LC and a steady state, and intermittency be- 
tween different LCs. A comparison between the experi- 
mental results and theory yields a partial agreement. 

Intermittency is a phenomenon in which a system re- 
sponse remains steady for periods of time (the laminar 
phase) which are interrupted by irregular spurts of rel- 
atively large amplitude dynamics (the turbulent phase). 
It arises in certain deterministic systems that are near 
a bifurcation in which a steady response is destabilized 
or destroyed This phenomenon also occurs in noisy 
systems in which the laminar response has a weak point 
in its local basin of attraction and is randomly bumped 
across the basin threshold, and then ultimately reinjected 
back to the laminar state, and the process repeats. This 
latter type of bursting behavior, which is relevant to the 
present system, is observed to occur in many other sys- 
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FIG. 1: (A) Measurement setup. (B) Schematic layout of the 
device. (C) Optical microscope image of the microbridge. 



tems, including Rayleigh-Bnard convection llOl . acoustic 
instabilities ^lliLturbulent boundary layers Il2ll. semicon- 
ducting lasers bhnking quantum dots [ij], sensory 
neurons [l5l| and cardiac tissues fl6l|. The presence and 
level of noise has a significant effect on all such systems, 
since perturbations affect the triggering of the system out 
of the laminar phase [13, [3 ■ The mean duration times 
of the laminar phase for a certain class of these systems 
scales in a manner that depends on the bifurcation pa- 
rameter and the noise level [l^ [13] ■ A special feature of 
the present system is that it exhibits a very sharp transi- 
tion between two types of operating states, namely, nor- 
mal conducting (NC) and superconducting (SC), which is 
modeled by equations with discontinuous characteristics. 
While the deterministic behavior of such non-smooth sys- 
tems (at least of low order) is generally well understood, 
including local and global bifurcations [2l| , the effects of 
noise in such systems has not been considered. 

The present experiments are performed using the setup 
depicted in Fig. [T^A) [12] ■ The resonator is stimulated 
with a monochromatic pump tone having an angular fre- 
quency ujp and power Pp . The power reflected off the res- 
onator is amplified at room temperature and measured by 
using both, a spectrum analyzer in the frequency domain 
and an oscilloscope, tracking the reflected power enve- 
lope, in the time domain. All measurements are carried 
out while the device is fully immersed in liquid Helium. 
A simplified circuit layout of the device is illustrated in 
Fig. HKB). The resonator is formed as a stripline ring 
made of Niobium Nitride (NbN) deposited on a Sapphire 
wafer [2^ . [23| , and having a characteristic impedance of 
50 ri. A feedline, which is weakly coupled to the res- 
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onator, is employed for delivering the input and output 
signals. A microbridge is monolithically integrated into 
the structure of the ring [23 |. 

The dynamics of our system can be captured by two 
coupled equations of motion Q. Consider a resonator 
driven by a weakly coupled feedline carrying an incident 
coherent tone 6'" = 6Q"e~"^p*, where b^^ is constant com- 
plex amplitude and tOp is the driving angular frequency. 
The mode amplitude inside the resonator can be written 
as i3e~"^p*, where B {t) is a complex amplitude, which is 
assumed to vary slowly on a time scale of l/wp. In this 
approximation, the equation of motion of B reads 

— ^[^{u;p-u;o)-l]B-^^,b"' + ^\ (1) 

where ujq is the angular resonance frequency and 7 = 
7i + 72 , where 71 is the coupling coefhcient between the 
resonator and the feedline and 72 is the damping rate 
of the mode. The term c'° represents an input Gaussian 
noise, whose time autocorrelation function is given by 
(c'"(t)c""(t')) = GloqS (t - t'), where the constant G can 
be expressed in terms of the effective noise temperature 
Tes as G = (27/ Wo) {kBTeff/hujo). The microbridge heat 
balance equation reads 

C^^2hu;o^2a\B\^-H{T-To), (2) 
at 

where T is the temperature of the microbridge, C is the 
thermal heat capacity, a is the portion of the heating 
power applied to the microbridge relative to the total 
power dissipated in the resonator {0 < a < 1), H is the 
heat transfer coefficient, and Tq = 4.2 K is the tempera- 
ture of the coolant. 

Coupling between Eqs. (P) and ^ originates by the 
dependence of the parameters of the driven mode ojq, 
71, 72 and a on the resistance and inductance of the 
microbridge, which in turn depend on its temperature. 
We assume the simplest case, where this dependence is 
a step function that occurs at the critical temperature 
Tc ~ 10 K of the superconductor, namely wq, 7i, 72 and 
a take the values wos, 7is, 72s and as respectively for the 
SC phase (T < Tc) of the microbridge and won, 7in, 72n 
and a,! respectively for the NC phase {T > Tc). 

Solutions of steady state response to a monochro- 
matic excitation are found by seeking stationary so- 
lutions to Eqs. |T]) and ^ for the noiseless case 
c'" — 0. The system may have, in general, up to 
two locally stable steady states, corresponding to the 
SC and NC phases of the microbridge. The stability 
of each of these phases depend on the corresponding 
steady states values Bg = 1^/2716'"/ [i {up — wqs) — 7s] 
and Bn = iy/2^b"^ / [i {ujp - ujon) - In] [see Eq. (P)]. 
A SC steady state exists only if |_Bs|^ < where 
Es — H {Tc — Tq) /2fiujosl2sCts, whereas a NC steady 
state exists only if |i3n| > £"„ where E^ — 

(Tc — Tq) /2fiwon72nan- Consequently, four stability 
zones can be identified in the plane of pump power 
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FIG. 2: Stability zones in the ujp — Pp plane: SC monos- 
table (MS(SC)), NC monostable (MS(NC)), bistable (biS), 
and astable (aS) (gray colored) zones. The region where a 
stable LC exists is marked with dashed line. The inset shows 
the three operating points (A, B and C) at which the measure- 
ments and theoretical analysis shown in Fig. |3]and|4]are done. 
The following parameters were used in the numerical simula- 
tion: a;os/27r = 3.49 GHz, 71s = 0.025aJos, 72s = OMOloos, 
ujoi^/LOQa — 1.017aJos, 7in = 0.2^LjQs and 7211 — 0.60ajos, 
C = 15.4fJ/K, H/C = 4.62a;os, = 700 K. loj, = 0.997i.Jos, 
Pp = -25.34 dBm (A), -25.78 dBm (B), -26.75 dBm (C). 



Pp cx l^o'^l^ - pump frequency ujp (see Fig. (2) 's']. Two are 
mono-stable (MS) zones (MS(SC) and MS(NC)), where 
either the SC or the NC phases is locally stable, respec- 
tively. Another is a bistable zone (BiS), where both 
phases are locally stable. The third is an astable zone 
(aS), where none of the phases are locally stable. 

The task of finding LC solutions of Eqs. (H]) and can 
be greatly simplified by exploiting the fact that typically 
7 ^ H/C in our devices, namely, the dynamics of the 
mode amplitude B [Eq. ([T])] can be considered as slow in 
comparison with the one of the temperature T [Eq. 
In this limit one finds by employing an adiabatic approx- 
imation Q that the temperature T remains close to the 
instantaneous value given by Ti = Tq + 2^^072 /H 
for most of the time except of relatively short time in- 
tervals (on the order of C /H) right after each switching 
event between the SC and NC phases. Consequently, as 
can be seen from the example trajectories shown in Fig. 
[3] (A-1), transitions from SC to NC phase occur near 
the circle \B\^ = Es, whereas transitions from NC to SC 
phase occur near the circle \B\ = E^. 

The important features of the system's dynamics can 
be captured by constructing a ID map 25]. Consider the 
case where £"„ < E^ and the amplitude B lies initially 
on the circle = namely B = y^^e^'^" where 
X e [0, 1]. Furthermore, assume that initially the system 
is in the SC phase, namely, T < Tc and consequently B 
is attracted towards the point Bg. The ID map D (x) is 
obtained by tracking the time evolution of the system for 
the noiseless case (c'" = 0) until the next time it returns 
to the circle \B\ = £'„ to a point B = y/E^e^^^^^^^ where 
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FIG. 4: Numerical [panels (A-3), (B-3) and (C-3)] vs. ex- 
perimental [panels (A-4), (B-4) and (C-4)] time traces for the 
three operating points A, B and C respectively. 



FIG. 3: (color online) Resonator's dynamics. Subplots A, B 
and C correspond to the three operating points A, B and C re- 
spectively which are marked in the inset of Fig. [5] In subplot 
(A) only a LC is locally stable, in subplot (B) intermittency 
between a LC and a SC steady state occurs, whereas only 
a SC steady state is locally stable in subplot (C). In panels 
(A-1), (B-1) and (C-1), which show the time evolution in B 
plane, a plus sign labels Bs and a cross sign labels Bn- These 
points are shown for reference and correspond to fixed points 
of the dynamics only when they exist in their respective do- 
mains, as defined in the text. Trajectories that return to the 
inner circle \Bf = En are colored in blue (dark), and tra- 
jectories that end at Bs are colored in yellow (gray). Panels 
(A-2), (B-2) and (C-2) show the corresponding ID maps. 



D (x) G [0,1]. In the adiabatic limit this can be done 
using Eq. ^ only [without explicitly referring to Eq. ^] 
since switching to the NC phase in this case occurs when 
the trajectory intersects with the circle = Eg. Note 
that in the aS zone of operation all points on the circle 
\Bf = En return back to it after a finite time. However, 
this is not necessarily the case in the other stability zones. 
Therefore, we restrict the definition of the ID map D (x) 
only for points on the circle =£"11 that eventually 
return to it. Other points will have a trajectory that ends 
at a steady state (NC or SC). 

Any fixed point of the ID map, namely a point for 
which D (xq) = xq, represents a LC of the system. The 
LC is locally stable provided that \dD/dx\^^^^ < 1 [2^ . 
The region in the ujp - Pp plane in which a locally stable 
LC solution exists was determined using the parameters 
of our device and it is marked with dashed line in Fig. [2l 

Figure [3] shows noiseless behvior of the resonator for 
the three operating points A, B and C, which lie near 
the border between the aS region and the MS(SC) one, 
and are marked in the inset of Fig. [2l Figure |4] shows a 
comparison of experimental data and numerical simula- 



tion for these operating points. The sample parameters 
used in the numerical simulations and are listed in the 
caption of Fig. [2l were determined using the same meth- 
ods detailed in Ref. 0. 

Subplot (A) shows the behavior at operating point A, 
which lies inside the aS zone. In panel (A-1) sample tra- 
jectories in B plane arc shown. The resultant ID map, 
which is plotted in panel (A-2), has a single fixed point 
corresponding to a single locally stable LC. The time evo- 
lution seen in panel (A-3) was obtained by numerically 
integrating the coupled stochastic equations of motion 
(HI) and ^ . The trace is then compared to experimental 
data taken from the same working point [panel (A-4)] . 

At operating point B [see Fig [3] and [4] subplot (B)] co- 
existence of a LC and a SC steady state occurs. The 
LC corresponds to the locally stable fixed point of the 
ID map seen in panel (B-2). On the other hand, all 
initial points on the circle = E^ that never return 
to it evolve towards the SC steady state B^. Numer- 
ical time evolution and experimental data shows noise- 
induced transitions between the two metastable solutions 
[panels (B-3) and (B-4) respectively] . At operating point 
C [see Fig |3] and [4] subplot (C)] the LC has been anni- 
hilated by a discontinuity-induced bifurcation t211 and 
consequently only steady state response is observed. 

To further study noise-induced transitions we fixed LUp 
and vary Pp starting from MS(SC) zone Pp = — 26.7dBm 
to the aS zone Pp = —25.6 dBm [vertical line in Fig. 
and took 1/i sec long time traces of Prcf [similar to those 
seen in Fig. |l](A-4), (B-4) and (C-4)]. The average life- 
time of both LC and SC steady state, namely, the average 
time the system is in one solution before making a transi- 
tion to the other one, were determined from these traces. 
This data, compared to numerical simulation prediction 
is shown in Fig. [5l 

In another experiment using a similar device we ob- 
serve intermittency of two different LCs (see Fig. (S]). 
Panel (A) shows spectrum analyzer measurement of the 
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FIG. 5: Experimental data of life time of LC (circles) 
and steady state (pluses) compared to numerical prediction 
(solid and dashed lines respectively). Mesurements above 
— 25.8dBm saturate to 1/isec as this is the maximal mesuer- 
ment time. 



reflected power as a function of the pump power Pp. 
Two distinct LCs having frequencies /i ~ 60 MHz and 
/2 ~ 80 MHz are observed. For Pp < —33.5 dBm only a 
LC at frequency /i is visible. In the range —33.55 dBm < 
Pp < —33.35 dBm both LCs are seen, whereas for high 
pump power Pp > —33.5 dBm only a LC at frequency /2 
is seen. Panels (B) and (C), demonstrates the behav- 
ior in the intermediate region, where both LCs are ob- 
served in frequency and time domain respectively for 
Pp = —33.5 dBm [marked by a dashed line in panel (A)]. 




In general, intermittency of two (or more) different LCs 
can be theoretically reproduced using our simple model. 
However, we were unable to numerically obtain this be- 
havior without significantly varying some of the system's 
experimental parameters. This discrepancy between ex- 
perimental and theoretical results suggest that a further 
theoretical study is needed in order to develop a more 
realistic description of the system. 
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FIG. 6: Experimental demonstration of intermittency be- 
tween two LCs. Panel (A) shows the reflected power P^ef 
as a function of the offset frequency A/ = {ui — L0p)/2n and 
the pump power Pp. Panel (B) shows a cross section of panel 
(A) which is indicated by a dashed line. The frequencies /i 
and /2 of the two LCs are marked. Panel (C) shows a time 
trace of the reflected power taken at the same value of Pp. In 
this experiment LUp/2n = 6.61 GHz 
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